Domain walls of ferroelectric BaTiC^ within the Ginzburg-Landau-Devonshire 
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Mechanically compatible and electrically neutral domain walls in tetragonal, orthorhombic and 
rhombohedral ferroelectric phases of BaTi03 are systematically investigated in the framework of the 
phenomenological Ginzburg-Landau-Devonshire (GLD) model with parameters of Ref. [Hlinka and 
Marton, Phys. Rev. 74, 104104 (2006)]. Polarization and strain profiles within domain walls are 
calculated numerically and within an approximation leading to the quasi-one-dimensional analytic 
solutions applied previously to the ferroelectric walls of the tetragonal phase [W. Cao and L.E. 
Cross, Phys. Rev. 44, 5 (1991)]. Domain wall thicknesses and energy densities are estimated for all 
mechanically compatible and electrically neutral domain wall species in the entire temperature range 
of ferroelectric phases. The model suggests that the lowest energy walls in the orthorhombic phase 
of BaTiOs are the 90-degree and 60-degree walls. In the rhombohedral phase, the lowest energy 
walls are the 71-degree and 109-degree walls. All these ferroelastic walls have thickness below 1 nm 
except for the 90-degree wall in the tetragonal phase and the 60-degree S-wall in the orthorhombic 
phase, for which the larger thickness of the order of 5nm was found. The antiparallel walls of the 
rhombohedral phase have the largest energy and thus they are unlikely to occur. The calculation 
indicates that the lowest energy structure of the 109-degree wall and few other domain walls in the 
orthorhombic and rhombohedral phases resemble Bloch-like walls known from magnetism. 

PACS numbers: 77.80.-e, 77.80.Dj, 77.84.Dy 
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I. INTRODUCTION 



Domain structure is an important ingredient in func- 
tionality of ferroelectric materials. Among others, it 
has impact on their nonlinear optical properties, dielec- 
tric permittivity and polarization switching phenomena. 
Since domain boundaries in ferroelectric perovskite mate- 
rials can simultaneously play the role of the ferroelectric 
and ferroelastic walls, such domain walls also strongly 
influence the electromechanical material properties: they 
facilitate switching of spontaneous polarization and spon- 
taneous deformation, thus giving rise to a large extrinsic 
contribution to e.g. piezoelectric constants, which makes 
ferroelectric materials extremely attractive for applica- 
tions. The domain structure also provides additional de- 
gree of freedom for tuning of material properties. In gen- 
eral, further development of domain engineering strate- 
gies requires deeper understanding of the physics of fer- 
roelectric domain wall itself. 

The Ginzburg-Landau-Devonshire (GLD) theory pro- 
vides a feasible tool for such a purpose. Landau- 
Devonshire model describes phase-transition properties 
of single-domain crystal using a limited number of pa- 
rameters, which are determined experimentally (or re- 
cently also using ab-initio methods). Introduction of 
Ginzburg gradient term to the free energy functional 
enables addressing nonhomogeneous multi-domain fer- 
roelectric state. The GLD model was previously used 
for computation of domain wall properties in ferroelec- 
tric materials (e.g. Refs.[ll-H) and in phase-field com- 
puter modeling of domain formation and evolution^— 
The GLD model can be regarded as a bridge model cov- 
ering length-scales inaccessible by ab-initio and micro- 



mechanical models. 

BaTiOs represents a typical ferroelectric material 
which undergoes a sequence of phase transitions from 
high-temperature paraelectric cubic PrnZm {Oh) to 
the ferroelectric tetragonal PAmm (C^), orthorhombic 
Amm2 {C^v) and rhombohedral Rim {C^ v ) phase. En- 
ergetically equivalent directions of spontaneous polar- 
ization vector, identifying possible ferroelectric domain 
states in a particular ferroelectric phase, are displayed in 
Fig.[U Domain boundaries separating two domain states 
are characterized by the rotational angle needed to match 
spontaneous polarizations on both sides of the bound- 
ary. For example, the boundary separating domains with 
mutually perpendicular spontaneous polarization is com- 
monly called 90° wall, while the one between antiparal- 
lel spontaneous polarization regions is called 180° wall. 
Other angles are possible in orthorhombic and rhombohe- 
dral phases of BaTiOs, where the spontaneous polariza- 
tion is oriented along the cubic face diagonals and body 
diagonals, respectively. 

The aim of this paper is to calculate basic characteris- 
tics of all electrically neutral and mechanically compat- 
ible domain walls in all ferroelectric phases of BaTiOs . 
For a better comparison of domain wall properties like 
their thickness or energy density, we employ an Ising- 
like approximation leading to previously proposed an- 
alytically solvable one-dimensional solutions^ The pa- 
per is organized as follows. In Section II. we give an 
overview of the different kinds of mechanically compati- 
ble domain walls in the three ferroelectric phases. It fol- 
lows from general theory about macroscopic mechanical 
compatibility of adjacent domain statesj 10 ' 11 The GLD 
parameters used for calculation of the domain wall prop- 
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FIG. 1: Directions of spontaneous polarization in the (a) 
tetragonal, (b) orthorhombic and (c) rhombohedral phase. 
Angles between polarization direction '0' and its symmetry 
equivalent ones are indicated. 



erties in BaTiOa are the same as in our preceding work, 4 
but for the sake of convenience, the definition and the 
GLD model and its parameters are resumed in Section 
III. Section IV. is devoted to the description of the com- 
putational scheme and approximations applied here to 
solve analytically the Euler-Lagrange equations. The 
main result of our study - systematic numerical evalu- 
ation of thicknesses, energies, polarization profiles and 
other properties for different domain walls, is presented 
in Section V. Sections VI. and VII. are devoted to the dis- 
cussion of validity of used approximations and the final 
conclusion, respectively. 



II. MECHANICALLY COMPATIBLE DOMAIN 
WALLS IN BARIUM TITANATE 

The energy-degeneracy of different directions of spon- 
taneous polarization leads to the appearance of ferroelec- 
tric domain structure. Individual domains are separated 
by domain walls, where the polarization changes from 
one state to another. Here, only planar domain walls 
are considered. Orientations of mechanically compatible 
domain walls are determined by the equation for mechan- 
ically compatible interfaces separating two domains with 
the strain tensors oo) and e^oo) : 



E 



>(oo) 



j(-Oo)]3 



= 



(1) 



Systematic analysis of this equation using symmetry ar- 
guments has been done e.g. in Refs. ll(ML2l In general, 
the number N of mechanically compatible domain walls 
separating two particular domain states can have only 
one of the three values: N = 0, N — 2 or TV = oo. 
In case of N — 2 there exist two mutually perpendic- 
ular domain walls. Each of them is either a crystallo- 
graphic (W/-type) wall or non-crystallographic (5-type) 
wall. Orientation of the FF/-wall is fixed by symmetry of 
the crystal, while orientation of the S-wall is determined 
by components of the strain tensor in adjacent domains 
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FIG. 2: Set of mechanically compatible and electrically neu- 
tral domain walls in the three ferroelectric phases of BaTiOa . 
In the case of 180° domain walls, where the orientation is 
not determined by symmetry, walls with the most important 
crystallographic orientations are displayed. 



(and its orientation can be therefore dependent on tem- 
perature). For N — oo there exists infinite number of 
wall orientations, some of them may be preferred ener- 
getically. 

Further, the electrically neutral domain walls will 
be considered^ It implies that the difference P(oo) — 
P(— oo) between the spontaneous polarizations in the ad- 
jacent domains is perpendicular to the unit vector s, nor- 
mal to the domain wall: 



(P(oo)-P(-oo))-s = 0. 



(2) 



We also define a unit vector r || (P(oo) — P(— oo)), which 
identifies the component of the spontaneous polarization 
which reverses when crossing the wall. Then the charge 
neutrality condition @ can be expressed as r • s = 0. 
Finally, let us introduce a third base vector t = r x s, 
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which complements the symmetry- adapted ortho normal 
coordinate system (r, s, t). 

BaTiOa symmetry allows a variety of domain wallsi 1 ^ 
Ferroelectric walls of BaTi03 can be divided in two 
groups - the non- ferroelastic 11 walls separating domains 
with antiparallel polarization (e mn (oo) — e m „(— oo) = 0, 
N = oo) and the ferroelastic walls with other than 
180° between polarization in the adjacent domain states 
(JV = 2). The r ■ s = condition implies that the neu- 
tral non-ferroelastic walls are parallel to the spontaneous 
polarization, and the neutral ferroelastic walls realize a 
"head-to-tail" junction. The set of plausible neutral and 
mechanically compatible domain wall types are schemat- 
ically shown in Fig.[2j Domain walls are labeled by a 
symbol composed of the letter specifying the ferroelectric 
phase (T, O or R staying for the tetragonal, orthorhombic 
or rhombohedral, resp.), number indicating the polariza- 
tion rotation angle (180, 120, 109, 90, 71 or 60 degrees) 
and, if needed, the orientation of the domain wall normal 
with respect to the parent pseudo-cubic reference struc- 
ture. 

Our choice of the base vectors r, s and of the sponta- 
neous polarization and strain components in the adjacent 
domain pairs for each domain wall type shown in Fig. [2] 
are summarized in TableQ] Base vectors coincide with 
special crystallographical directions, except for the O60 
wall where the s and t vectors depend on the orthorhom- 
bic spontaneous strain (see Tabled as follows: 5 
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with L>i = V2D 2 , D 2 = ^(e a - e c ) 2 + 2eg and with 
e a , eb, and e c defined in Table I. 

Although only the neutral walls are discussed in the 
following, the Fig. [5] is actually helpful in enumeration of 
all possible mechanically compatible domain wall species 
in BaTiOa. In principle, mechanical compatibility allows 
180° VFoo-type domain walls with an arbitrary orienta- 
tion of the domain wall in all three ferroelectric phases 
(T180, O180, R180). Obviously, they are electrically neu- 
tral only if the domain wall normal is parallel with the 
spontaneous polarization. Ferroelastic walls exists in mu- 
tually perpendicular pairs. In the tetragonal phase, there 
exist 90° Wf -type domain walls (T90), either charged 
(head-to-head or tail-to-tail) or neutral (head-to-tail). 
The orthorhombic phase is more complex. In the case 
of 60° angle between polarization directions, the N=2 
pair is formed by a charged W/-type wall and neutral S- 
type wall. The case of 120° angle is similar but W/-wall 
is neutral and S'-wall is charged. In addition, there are 
again charged or neutral 90° Wy-walls (O90). The rhom- 
bohedral phase has pairs of charged and neutral Wf-type 
domain walls with the angle between polarizations either 



109° or 71° (R109 or R71, resp.). Since only neutral 
walls are discussed here, the 5-type domain wall will be 
referred to as O60 and Wf wall as O120. 



III. GLD MODEL FOR BARIUM TITANATE 

Calculations presented in this paper are based on the 
GLD model with anisotropic gradient terms, reviewed in 
Ref.0- The free energy F is expressed in terms of polar- 
ization and strain field taken for primary and secondary 
order-parameter, resp.: 



F[{P i ,P iij ,e ij }]= / /dr, 



(4) 



where the free energy density / consists of Landau, gra- 
dient, elastic and electrostriction part 



I = fi e) {P l } 



fc{P, e«} + f q {P, e tJ } + f G {P id }. (5) 



The Landau potential considered here is expanded up 
to the sixth order in components of polarization for the 
cubic symmetry (0\): 
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with the three temperature-dependent coefficients 
ai, an and am, as in Ref.fill This expansion produces 
the 6 equivalent domain states in the tetragonal phase, 
12 in the orthorhombic, and 8 in the rhombohedral phase 
(see Fig.HJ). 

Dependence of the free energy on the strain is en- 
countered by including elastic and linear-quadratic elec- 
trostriction functionals Fq and F q , resp. Their corre- 
sponding free energy densities are 



fc 



2 ^pCpo-Ga 



and 



fq — —qijkiSijPkPi 
i (dui/dxj + duj/dxi) and C a p, q a p 



(7) 



(8) 



where 

components of elastic and electrostriction tensor in Voigt 
notation, Cu = Cmij C\ 2 = C1122, C44 = C1212, <Zn = 
<7mi, q\ 2 = <7ii22, but 1744 = 2<7n22- 

The elastic and electrostriction terms result in re- 
normalization of the bar expansion coefficients ot[} and 
when minimizing the free energy with respect to 
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TABLE I: Cartesian components of switching vectors r, domain wall normals s, and boundary conditions for polarization and 
strain in adjacent domain states for the inspected domain walls. Vector S060 is defined in Eqn. @, Po stands for magnitude 
of spontaneous polarization. As usual, spontaneous quantities are those minimizing GLD functional. Numerical values used in 
this work are given in Sec'IVl 

Wall r s P(-oo)/P P(oo)/P e(-oo) e(ooj 

T180{001} (1,0,0) (0,0,1) (1,0,0) (-1,0,0) (ejj , ex, el, 0, 0, 0) (ejj , e ± , e±, 0, 0, 0) 

T180{011} (1,0,0) (0,^5,7=) (1,0,0) (-1,0,0) (en.ej., ex, 0,0,0) (ey, ex, ex, 0,0,0) 

T90 (75'75' ) (75'^' ) t 1 ' ' ) (0'- 1 ' ) (e,|, e x , ex, 0,0,0) (ex, e,| , e x , 0, 0, 0) 

O180{110} (^5,^5,0) (7^,75,0) (7^,7^,0) (75, 7^0) (e a ,e a ,e c ,0,0,2 eb ) (e a , e a , e c , 0, 0, 2e b ) 

O180{001} (75,75,0) (0,0,1) (75,75,0) (75,75,0) (ea,e a , ec ,0,0,2e b ) (e a , e a , e c , 0, 0, 2e b ) 

O90 (0,1,0) (1,0,0) (75,75,0) (75,75,0) (e a , ea ,ec,0,0,2 eb ) (e a , e a , e c , 0, 0, -2e b ) 

060 (73' '72) 8060 (72' 75'°) (°'72'7f) (ea,e a ,e c ,0,0,2e b ) (e c , e a , e a , -2e b , 0, 0) 

0120 (Tif'Te'vf) (75'°'7l) (75'75'°) ( >v3'?f) (ea,e a ,e c ,0,0,2e b ) (e c , e a , e a , -2e b , 0, 0) 

R180{ll0} (7^,7^,7^) (^5,75,0) (75,7^,7^) (tI.^.tI) (e a ,e a) e a ,2e b ,2e b ,2e b ) (e a , e a , e a , 2e b , 2e b , 2e b ) 

R180{211} (7^7=) (§,7^) (73.7f.73) ($,4,$) (e a , ea , Ca ,2 eb ,2 eb ,2e b ) (e a , e a , e a , 2e b , 2e b , 2e b ) 

R109 (75'73' ) (O' ' 1 ) (TI'Tf'Tf) (TI'TI'Tf) (ea, e a , e a , 2e b , 2e b , 2e b ) (e a , e a , e a , -2e b , -2e b , 2e b ) 

R71 (0,1,0) (75,0,75) (71.75:7?) (75.7g.7g) (e a ,e a ,e a ,2e b ,2e b ,2e b ) (e a , e a , e a , -2e b , 2e b , -2e b ) 



0) 



strains (in the homogeneous sample). The bar a n , a 12 
and the relaxed an, a±2 coefficients are related asi^ 
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The Ginzburg gradient term /g is considered in the 
form 
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It was pointed oul4 that the tensor of gradient constants 
of BaTiOa is highly anisotropic with fundamental conse- 
quences on predicted domain wall properties. Up to now, 
the isotropic gradient tensor Gij was mostly employed 
in the computations. 

The material-specific coefficients in the model are as- 
sumed being constant, except for the three Landau po- 
tential coefficients 

a t = 3.34 x 10 5 (T- 381) 
an = 4.69 x 10 6 (T- 393) - 2.02 x 10 8 
a m = -5.52 x 10 7 (T- 393) + 2.76 x 10 9 , (12) 



where T is absolute temperature^ The phase transi- 
tions occur in this model at the temperatures 392.3 K 
(C->T), 282.5K (T->0), and 201.8K (O^R)^ All 
phase transitions are of the first order, the local min- 
ima corresponding to the tetragonal, orthorhombic and 
rhombohedral phase exist for this Landau potential 
between 237 K and 393 K, between 104 K and 303 K, 
and below 256 K, respectively. Full set of temperature 
independent parameters of the GLD model reads^^ 
a i2 = 3.230 x 10* Jm 5 C" 4 , a n2 = 4.470 x 10 9 Jm 9 C- 6 
"123 = 4.9 1 x 10 9 Jm 9 C- 6 , Gn = 51 x 10~ n Jm 3 C~ 2 
G12 = -2 x 10- 11 Jm 3 C- 2 , G44 = 2 x 10- 11 Jm 3 C- 2 
gn = 14.20 x 10 9 JmC- 2 , 912 = -0.74 x 10 9 JmC- 2 
944 = 1.57 x 10 9 JmC- 2 , Gn = 27.50 x 10 10 J m - 3 
G12 = 17.90 x 10 10 Jm- 3 and G44 = 5.43 x 10 10 JiTr 3 . 



IV. STRAIGHT POLARIZATION PATH 
APPROXIMATION 

Let us consider a single mechanically compatible and 
electrically neutral domain wall in a perfect infinite 
stress-free crystal. Within the continuum GLD theory, 
such domain wall is associated with a planar kink solution 
of the Euler-Lagrange equations of the GLD functional 
for polarization vector and the strain tensor.i 2 -^^ Do- 
main wall type is specified by selection of the wall normal 
s and by the two domain states at s = —00 and +00. This 
ideal geometry implies that the polarization and strain 
vary only along normal to the wall s and domain wall 
can be thus considered as a trajectory in order-parameter 
space. 

Even if the domain wall is neutral, in its central 
part the local electric charges can occur due to the 
position-dependent polarization. The additional assump- 
tion V • P = ensures the absence of charges in the 
whole domain wall. Then the local electric field is zero 
and the electrostatic contribution vanishes. Such strictly 
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charge-free solutions^ were found to be an excellent ap- 
proximation for ideal dielectric materials^ The condition 
V • P = implies that the polarization vector variation 
is restricted to a plane perpendicular to s (the trajectory 
is constrained to P s = P s (±oo) plane). 

In fact, the polarization trajectories representing the 
T180 and T90 walls calculated under VP = constraint 
were found to be the straight lines connecting the bound- 
ary values This greatly simplifies the algebra and 
the solution of the variational problem can be found an- 
alytically. Therefore, we decided to impose the condition 
of a direct, straight polarization trajectory for the varia- 
tional problem of all domain wall species of BaTiC>3 . This 
condition, further referred as straight polarization path 
(SPP) approximation, implies that both P s and Pt po- 
larization components are constant across the wall. We 
shall come back to the meaning and possible drawbacks 
of this approximation in Section VI. 

Polarization and strain in the mechanically compat- 
ible and electrically neutral SPP walls can be cast 
in the form: P = [P T (s), P s (±oo), P t (±oo)] and e = 
[e rr (±oo), e ss (s), e tt (±oo), 2e st (s), 2e rt (±oo), 2e rs (s)], re- 
spectively. The s-dependent strain components are cal- 
culated from the mechanical equilibrium condition 
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4-f dxj 

3=1 J 
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dxj V de 
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(13) 



Boundary conditions for stress and the fact that all quan- 
tities vary only along direction s imply 



den 



C = 



(14) 



for ij G {ss, st, rs} with the integration constant C to be 
determined from boundary conditions. 

The Euler-Lagrange equation for polarization reduces 

to 



d Of 8f 



ds dP T . s dP T 



= . 



(15) 



where the strain components from Eqn. (|14[) were substi- 
tuted into /. Thus, the elastic field was eliminated. It 
turns out that the resulting Euler-Lagrange equation for 
P T is possible to rewrite in the form 



9- 



d 2 p(s) 
' ds 2 



2a lP (s) + 4a llP 3 (s) + 6a m p 5 (s) , (16) 



where p(s) stands for P r (s), and where the coefficients 
g, ai, an and am, different for each domain type, depend 
only on the material tensors. The boundary values are 

-Pr(oo) = Pr(-Oo) = p^. 

Solution of the Euler-Lagrange equation (Til))) is well 
known»iii£~— We shall follow the procedure of Ref.l20l. 
Integrating Eqn. (ITfJl) one can obtain the equation: 



f(|) =f*L<P) 
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FIG. 3: Profile of the "reversed" polarization component 
P r (Eqn. (23) for a domain wall with the shape factor A = 1 
(solid) and A — 3 (dashed). Both profiles have the same 
derivative for s/2£ = (dotted) and therefore also the same 
thickness (indicated by vertical lines) according to the defini- 
tion in Eqn. ([23)). 



where (see Ref.|j) 

f EL (p) = aip(s) 2 + anp(s) 4 + ainp(s) 6 . (18) 

The function /el(p) is a double- well "Euler-Lagrange" 
potential with two minima ±Poo , where 



P~ 



-an + yafi - 3flmai 
3am 



(19) 



The differential Eqn. (116)) has the analytical solution 

smh(s/£') 



P(s) =Po 



^L + sinh 2 (s/£') ' 



(20) 



where 



and 



3«lllPoo + Qn 

2amP 2 c + an 



(21) 



(22) 



Quantity A determines deviation of the profile (|2T)|) from 
the tanh profile, which occurs for the 4" l -order potential 
(i.e., am = 0, A = 1). 

The domain wall thickness (Fig. El is defined as 



(23) 



with U being the energy barrier between the domain 
states: 




U = /bl(0) - /bl(poo) = 2011^ + anp^. 



(24) 
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The surface energy density of the domain wall is 

E = / Uel{s) - fEL(Poo))ds 



- Poo ^2gU[A 5 / 2 I(A) 



where 



cosh 2 (h)dh 



(25) 



(26) 



(A + cosh 2 (/i) - l) 3 ' 

The domain wall characteristics depend on the coeffi- 
cients g, ax, an, am through Eqs. §M [201 121 dSJ) ■ The 
explicit expressions for these coefficients are summarized 
for various domain walls in Table II. The expressions are 
simplified using the notation inspired by Ref.Q]. For all 
phases we are using: 
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for tetragonal and orthorhombic phases we are abbrevi- 
ating 
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(27) 



The expressions for the coefficients of the T180 and 
T90 domain walls are equivalent to the previously pub- 
lished expressions Derivations for O60, O120, and 
RI80{2I1} walls lead to complicated formulas, and there- 
fore only numerical results for g,a\,a\\, and am coeffi- 
cients are presented here. 
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FIG. 4: Dependence of correction factor in the expression for 
domain wall energy density (25) as a function of the shape 
coefficient A. Full points indicate values of A for particular 
domain walls considered in Table III. 
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V. QUANTITATIVE RESULTS 

Advantage of the GLD approach is that domain wall 
properties can be obtained at any temperature. In Ta- 
ble lHIl we give numerical results for domain wall param- 
eters at one particular temperature for each of the ferro- 
electric phases: at 298 K, 208 K, and 118 K for the tetrag- 
onal, orthorhombic, and rhombohedral phase, resp. Cor- 
responding numerical values of the spontaneous quanti- 
ties appearing in TableQ] are: Pq = 0.265 Cm -2 , en = 
Q u p 2 = 7.77 x 10~ 3 and e±_ = Qi 2 P 2 = -3.18 x 10~ 3 
(the tetragonal phase, 298 K); P = 0.331 Cm" 2 , e a = 
±Qi2 P 2 _nsv m-3 ; ec = Q X2 p2 = -4.96 x 10~ 3 

0.79 x 10 -3 (the orthorhombic 
Pa = 0.381 Cm" 2 , e a = Qll+2Q 



, P 2 = 3.58 x 10" 

_ Q44 p2 



and eb 
phase, 208 K 
0.97 x 10~ 3 and e b 



Qi 



p2 _ 



■P 2 



0.70 x 10~ 3 (the rhom- 
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TABLE II: Parameters characterizing various mechanically compatible and neutral domain wall species of BaTi03-like fer- 
roelectrics within the SPP treatment described in the Section IV. Results for O60, O120, and R180{211} walls as well as a 
few other parameters are omitted because the corresponding analytical expressions through the GLD model parameters and 
spontaneous values of polarization and strain are too complicated. 
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TABLE III: Predicted values of thickness and planar energy density of domain wall species illustrated in Fig. 2 together with the 
determining parameters appearing in the SPP treatment. Results are evaluated from the BaTiOs-specific GLD model at three 
selected temperatures corresponding to the tetragonal, orthorhombic and rhombohedral phase, respectively. Domain walls for 
which relaxing of the Pt component results in a lower-energy CPP solution (see in Section VI.) are denoted by Numerical 
values are in SI units (2£ in nm; S in mj/m 2 ; U in MJ/m 3 ; poo in C/m 2 ; g in 10 -11 kgm J s~ 2 C~ 2 ; ai in 10 7 kgm 3 s~ 2 C~ 2 ; 
an in 10 8 kgm 7 s- 2 C" 4 ; am in 10 9 kgm n s- 2 C" 6 ). 



Domain wall 




£ 


U 


A 


Poo 


9 


01 


an 


am 








Tetraj 


jonal phase 


(298 K) 










T180{001} 


0.63 


5.9 


6.41 


1.43 


0.265 


2.0 


-14.26 


1.69 


8.00 


T180{011} 


0.63 


5.9 


6.41 


1.43 


0.265 


2.0 


-14.26 


1.69 


8.00 


T90 


3.59 


7.0 


1.45 


1.09 


0.188 


26.5 


-7.86 


9.53 


3.12 








Orthorhombic phase (208 K) 










oisoiiiojt 


2.66 


31.0 


8.20 


1.70 


0.331 


26.5 


-9.71 


-2.75 


4.36 


O180{001} 


0.70 


8.9 


8.95 


1.64 


0.331 


2.0 


-11.07 


-2.13 


4.36 


O90 


0.72 


4.3 


4.26 


1.50 


0.234 


2.0 


-11.62 


-0.08 


12.97 


O60 


3.62 


5.3 


1.09 


1.08 


0.166 


26.1 


-7.64 


12.14 


4.36 


O120 1 " 


1.70 


13.7 


5.82 


1.47 


0.287 


10.2 


-10.82 


0.50 


4.92 








Rhombohedral phase (118 K) 










RISO^TO} 1 ^ 


2.13 


36.0 


11.81 


1.83 


0.381 


18.3 


-9.53 


-3.64 


3.17 


R180{211} t 


2.13 


36.0 


11.81 


1.83 


0.381 


18.3 


-9.53 


-3.64 


3.17 


R109+ 


0.70 


7.8 


7.81 


1.65 


0.311 


2.0 


-10.84 


-2.56 


5.60 


R71 


0.74 


3.7 


3.52 


1.58 


0.220 


2.0 


-10.31 


-2.42 


17.94 



8 




100 200 300 400 100 200 300 400 

T(K) T(K) 

FIG. 6: Temperature dependence of the thickness (a) and energy density (b) for various mechanically compatible and neutral 
domain walls of BaTiC>3 estimated within SPP approximation. The SPP values for domain walls for which relaxing of the Pt 
component results in a lower-energy CPP solution (see in Section VI.) are shown by dashed lines. Thickness of O180{001}, 
O90 has a very close temperature dependence. The results for T180 as well as for R180 walls are not distinguished by specific 
orientations of the wall normal, since for both cases the angular dependence of thickness and energy of SPP solutions on the 
domain wall normal is negligible (see Tab. lIII|) . Vertical lines mark phase-transition temperatures. 



bohcdral phase, 118 K). The base vectors for the O60 
domain wall (see Eqn.[3]) at 208 K are 

r O60 = (0.707,0,0.707) 

s O60 = (0.701,0.130,-0.701) 

toeo = (-0.092,0.991,0.092) . (28) 

The right four columns of the Table lLUl contain corre- 
sponding numerical values of the domain wall coefficients 
g, di, an, and am derived from GLD parameters and 
spontaneous order-parameter values using above derived 
analytical expressions, mostly given explicitly in TablelLTl 
The left part of the Table lLUl contains the key domain wall 
properties such as wall thickness 2£ and energy density 
E. 

Clearly, T90 and O60 domain walls are considerably 
broader then others. T90 wall in the tetragonal phase is 
predicted to be 3.59 nm thick at 298 K (as already cal- 
culated in Ref.Hj) and S-type O60 domain wall in the 
orthorhombic phase is predicted to be 3.62 nm thick at 
208 K. Since the wall thickness is greater than the lat- 
tice spacing, the pinning of the walls is weak^i and they 



can be easily moved. Moreover, in case of O60 wall, the 
pinning is further suppressed due to 'incommensurate' 
character of Miller indices of the wall normal. 

As follows from Eqn. (|2"3"1) . domain wall thickness is 
determined by quantities g, U, and Poo. By inspec- 
tion of their values in Table HIIl we see that the most 
important factor is the coefficient g. Indeed, the do- 
main walls R71, R109, O90, O180{001}, T180 with 
g = 2x 10~ n kgm 5 s -2 C -2 are all very narrow (thickness 
below 1 nm) , and the thickness of various walls mono- 
tonically increases with increasing value of g. Let us 
stress that coefficients di, an, am, and g depend on 
the direction of the domain wall normal. For exam- 
ple, 0180(110} with g = 26.5 x 10" 11 kgm 5 s- 2 CT 2 
is almost four times broader than O180{001} with g = 
2 x 10- 11 kgm 5 s- 2 C- 2 . 

In the absence of other constraints, the probability of 
appearance of domain wall species should be determined 
by the surface energy density E. Therefore, in the or- 
thorhombic phase, the thinner O180{001} wall is more 
likely to occur than the 0180(110} one. Interestingly, 
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the O180{001} wall has almost the same thickness as the 
O90 wall, while in the tetragonal phase, it is the 90° wall 
which is much thicker than 180° wall (3.59 nm compared 
to 0.63 nm). 

In general, the normal of the neutral 180°-domain wall 
can take any direction perpendicular to the spontaneous 
polarization of the adjacent domains. Therefore also E 
depends on the orientation of the wall normal (in the s — t 
plane). We have checked in the orthorhombic phase that 
the 0180{001} and 0180{110} correspond to the extremes 
in the angular dependence of E, which is monotonous 
between them. This strong angular dependence is corre- 
lated with the anisotropy of the tensor Gijki- However, 
g is independent of the wall direction in the tetragonal 
and rhombohedral phase, and also the variation of coef- 
ficients Oi, On, an d am is insignificantly small, so that 
the effective directional dependence of 180° domain wall 
properties in these phases is negligible. 



T1 80(001} 



T180{011} 



T90 



The values of the shape factor A (see Table llllj) deter- 
mines the deviation of the P r (s) polarization profile of the 
domain walls from the simple tanh form. For all studied 
cases the values of A range between 1 and 2, where the 
correction factor A 5 / 2 I(A) appearing in the Eqn. (|25l) is 
almost linear function of A, as it can be seen in Fig. 2) It 
means that the shape deviations are much smaller than 
those shown by broken line in Fig.[3l Unfortunately, in 
the case of broad walls T90 and O60, which are good 
candidates for study of the structure of the wall central 
part, A is almost one. 

Eqn. (TH)) can be also used to evaluate local strain vari- 
ation in the wall. In Fig. [5] the profiles of strain tensor 
components for T90 domain wall are shown for illustra- 
tion. Obviously, 633,623, and e±3 strain components are 
strictly constant and equal to their boundary values as it 
follows from mechanical compatibility conditions. The 
en and e22 components and polarization P r vary be- 
tween their spontaneous values. Let us stress that the 
're-entrant' shear component ei2 (it has the same value 
in both adjacent domains) approaches the non-zero value 
of about 5xl0~ 4 in the middle of the domain wall. 
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FIG. 7: Equipotential contours of the Euler-Lagrange poten- 
tials showing ELPS associated with the P S =P S (±00) plane 
for a set of BaTiOs domain walls considered in Table III. Po- 
larization scale given for the bottom left inset (in Cm" 2 ) is 
equally valid for all shown ELPS's. Energetically most favor- 
able domain wall solutions with trajectories restricted to the 
P s =P s (±oo) plane are indicated with bold lines. The solu- 
tions corresponding to a linear segment are denoted in the 
text as SPP walls, as opposed to the CPP solutions (see in 
Section VI.), which have a curved order-parameter trajectory. 



The temperature dependence of the thickness and sur- 
face energy density of the 12 studied domain wall species 
is plotted in Fig.lBh, and Fig.|Bb, respectively. Although 
some properties do vary considerably, e.g. thickness of 
T90 wall in the vicinity of the paraelectric-ferroelectric 
phase transition, the sequence of the thickness values as 
well as the surface energy-density values of different do- 
main wall types remain conserved within each phase. 

The temperature dependence of domain wall thick- 
ness follows the trend given by Eqn. (|2"3")l . It increases 
with increasing temperature due to the dependence on 
Poo (see dependence of U on in Eqn. PHO . Such 
behavior is well known also from the experimental 
observations. 3,22,23 



VI. CURVED POLARIZATION PATH 
SOLUTIONS 

So far we have investigated domain walls within the 
SPP approximation, i.e. components P s and P t , which 
are the same in both domain states, were kept constant 
inside the whole wall. This is quite usual assumption 
made for a ferroelectric domain wall. Nevertheless, the 
full variational problem, where all three components of 
P could vary along s coordinate, leads in general to a 
lower energy solution corresponding to a curved polar- 
ization path (CPP) in the 3D primary order-parameter 
space. Appearance of the non-zero 're-entrant' compo- 
nents within a ferroelectric domain wall was considered 
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e.g. in works of Refe-HHH! For 180° walls, the polar- 
ization profile associated with SPP is sometimes denoted 
as the Ising-type wall. In contrast, the CPP solutions 
with non-zero P s , Pt are often considered as Neel and 
Bloch-like^ even though, in contrast with magnetism, 
the modulus of P is far from being conserved along the 
wall normal s. 

We have previously considered non-constant P s com- 
ponent of polarization in the 90° wall with explicit treat- 
ment of electrostatic interaction and realized that the 
deviations from SPP approximation are quite negligible. 4 
In general, non-constant P s would lead to non-vanishing 
V • P and finite local charge density, which in a perfect 
dielectric causes a severe energy penalty. The same sit- 
uation is expected for all domain wall species. However, 
there is no such penalty for non-constant P t -solutions. It 
was previously argued^ that the Bloch-like (with a con- 
siderable magnitude of Pt at the domain wall center) so- 
lutions could occur in orthorhombic BaTi03. Therefore, 
it is interesting to systematically check for existence of 
such solutions using our model. 

In order to study such Bloch-like solutions, we have cal- 
culated Euler-Lagrange potential in the order-parameter 
plane P s = P s (±oo) by integrating Euler-Lagrange equa- 
tions (Eqn.rni] and [H)) for all domain wall species from 
Table lllll similarly as e.g. in the Refs.Q]]!. Resulting 2D 
Euler-Lagrange potential surfaces (ELPS) are displayed 
in Fig. [3 In each ELPS, the bold lines indicate numeri- 
cally obtained domain-wall solution with the lowest en- 
ergy. The spatial step was chosen as 0.1 nm, P r and Pt 
were fixed to boundary conditions in sufficient distance 
from domain wall (6 nm) and initial conditions for P t 
were chosen so that the polarization path bypasses the 
energy maximum of the ELPS, and the system was re- 
laxed to the equilibrium. 

Among the twelve treated wall species, there are six 
cases where only the SPP solutions with P t =const exist: 
T180{001}, T180{011}, T90, O180{110}, O90, and R71. 
These solutions are clearly "Ising-like" . In all these cases, 
the (P r = 0, Pt = 0) point is the only saddle point of the 
ELPS. In the other six cases - O180{001}, O60, O120, 
R180{110}, R180{211}, and R109 - the ELPS has a max- 
imum at the (P r = 0, Pt =0), and the lowest energy solu- 
tions correspond to curved polarization paths. This sug- 
gests that the previously discussed SPP description may 
not necessarily be the proper approximation for these 
walls. Nevertheless, in the case of O60 and O120 walls 
the deviations from the SPP model are marginal, and 
only the remaining four solutions exhibit strong Bloch- 
like behavior. Moreover, the energy differences between 
SPP and CPP solutions were found to be almost negligi- 
ble, except for R180, where the CPP energy is by about 
10% lower in the entire temperature range of stability of 
the rhombohedral phase. Therefore, it is quite possible 
that in the case of O180{001}, R180{110}, R180{211}, 
and R109 walls both Bloch-like and Ising-like solutions 
may be realized. 

The deviations from SPP in the case of the almost 
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FIG. 8: Predicted profiles of polarization components for O60 
domain wall in BaTiOs at T = 208 K. Full line stands for the 
P r component of polarization vector, broken line for the Pt 
one. Analytically obtained SPP solution (a) practically does 
not differ from the numerically obtained CPP solution (b) 
with an unconstrained Pt component in this case. 



Ising-like O60 and O120 walls are associated with the 
fact that ELPS is not symmetric with respect to Pt = 
P t (±oo) mirror plane. In these cases not only the po- 
larization path and wall energies, but also domain wall 
thicknesses of the SPP and CPP counterpart solutions 
are very similar. This is demonstrated in Fig. [SI which 
shows polarization profiles of both SPP and CPP solu- 
tions for the O60 wall. 

Much more pronounced difference between domain 
wall profiles of SPP and CPP solutions are found in case 
of O180{001}, R180{110}, R180{211}, and R109 Bloch- 
like walls where the CPP trajectories bypass the (P r = 0, 
Pt = 0) maximum near the additional minima, which 
originate from 'intermediate' domain states either of the 
same phase or even of the different ferroelectric phases. 
In these cases, the inadequacy of SPP approximation is 
obvious. For example, the CPP of R180{211} wall seems 
to pass through a additional minimum corresponding to 
an 'orthorhombic' polarization state (consult correspond- 
ing inset in the Fig. [2]). As expected, the profile of such 
CPP solution deviates strongly from tanh shape and even 
definition of the wall thickness would be problematic (see 
Fig.©. 
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FIG. 9: Profiles of polarization components in R180{110} 
domain wall showing the SPP stationary trajectory (a) as well 
as the lower energy Bloch-like solution (b), both for the model 
parameters corresponding to T = 118 K. Full line stands for 
the P r component of polarization vector, broken line for the 
Pt one. The bottom panel demonstrates that Bloch solution 
may considerable modify the domain domain wall profile. 

VII. CONCLUSION 

The work reports detailed study of mechanically com- 
patible and electrically neutral domain walls in BaTiC>3. 
The investigation was done within the framework of the 
GLD model. Using the SPP approximation it was possi- 
ble to compare properties of various kinds of domain wall 
species from the same perspective. 

The phenomenological nature of the GLD model al- 
lowed to predict the temperature dependence of the do- 
main wall characteristics in the whole temperature range 
of ferroelectric phases. Its continuous nature gave us 
even the opportunity to deal conveniently with the non- 
crystallographic S-type domain wall, which has a general 
orientation with respect to the crystal lattice, and which 
is therefore difficult to cope with in discrete models rely- 
ing on periodic boundary conditions. 

The S-wall in the orthorhombic, as well as the 90° 
wall in tetragonal phase were both found to be about 
4nm thick and consequently are expected to be mobile, 



i.e. they could be easily driven by external fields, and 
they may thus significantly contribute to the dielectric 
or piezoelectric response of the material. 

For several temperatures, we have numerically investi- 
gated domain walls allowing for more complicated CPP 
solutions with non-constant Pf We have identified solu- 
tions, which could be considered as analogues of Bloch 
walls known from magnetism. Interestingly, in contrast 
with Ref.i, our model predicts the Ising-type profile of 
the O180{110} wall. At the same time the Bloch-like 
structure of the O180{001} wall is predicted. 

We believe that this kind of somewhat exotic walls ac- 
tually represent important generic examples of ferroelec- 
tric domain species, which should be anticipated in all 
ferroelectrics with several equivalent domain states dis- 
tinguished simultaneously by the orientation of the spon- 
taneous polarization and strain. They should be certainly 
taken into account in investigations of domain-wall phe- 
nomena in ferroelectric perovskites. At the same time, 
the energy differences between the Bloch-like and Ising- 
like solutions are rather subtle here. Therefore, in spite 
of the fairly good agreement for 180° and 90° domain 
walls between ab-initio calculations and predictions of 
this model in the tetragonal phased the preference for 
the calculated Bloch-like trajectories may not necessar- 
ily reproduced for the domain walls encountered in real 
BaTiC>3 crystal, since there is obviously a considerable 
uncertainty in the adopted material-specific GLD param- 
eters. In addition, predictions for the domain walls with 
very small thickness must be considered with a particu- 
lar caution since the description of the sharp domain wall 
profiles obviously touches the limits of the applicability 
of the continuous model. 

In conclusion, we have derived a number of qualita- 
tive and quantitative predictions for mechanically com- 
patible neutral domain walls of tetragonal, orthorhombic 
and rhombohedral BaTi03 on the basis of the previously 
proposed material-specific GLD model. We believe that 
the insight into the domain wall properties mediated by 
the provided analytical and numerical analysis could be 
helpful for understanding of domain wall phenomena in 
BaTiOa as well as in some other intensively investigated 
members of the ferroelectric perovskite family with same 
sort of macroscopic ferroelectric phases, for example in 
KNb0 3 , BiFe0 3 , PbTi0 3 or even PZT and perovskite 
relaxor-related materials. 
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